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Abstract 

We study the two-dimensional Gross-Pitaevskii theory of a rotating Bose gas in a disc-shaped trap 
with Dirichlet boundary conditions, generahzing and extending previous resuhs that were obtained 
under Neumann boundary conditions. The focus is on the energy asymptotics, vorticity and quah- 
tative properties of the minimizers in the parameter range | loge| Q, < e~^| Ioge|~^ where Q is the 
rotational velocity and the coupling parameter is written as with e ^ 1. Three critical speeds 
can be identified. At f2 = ~ I log £ I vortices start to appear and for |loge| ^ n < f2c2 ~ 
the vorticity is uniformly distributed over the disc. For SI > the centrifugal forces create a hole 
around the center with strongly depleted density. For Q <^ e^^|log£|^^ vorticity is still uniformly 
distributed in an annulus containing the bulk of the density, but at Q, = Q,C3 ~ e~^| loge|^^ there 
is a transition to a giant vortex state where the vorticity disappears from the bulk. The energy is 
then well approximated by a trial function that is an eigenfunction of angular momentum but one 
of our results is that the true minimizers break rotational symmetry in the whole parameter range, 
including the giant vortex phase. 

MSG: 35Q55,47J30,76M23. PACS: 03.75.Hh, 47.32.-y, 47.37.-(-q. 
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1 Introduction and Main Results 

The Gross-Pitacvskii (GP) theory is the most commonly used model to describe the behavior of rotating 
superfluids. Since the nucleation of quantized vortices is a signature of the superfluid behavior it is of 
great interest to understand that phenomenon in the framework of the GP theory. A fascinating example 
of superfluid is provided by a cold Bose gas forming a Bose-Einstein condensate (BEG). The possibility 
to nucleate quantized vortices in a rotating BEG has triggered a lot of interest in the last decade, both 
experimental and theoretical (see the reviews |Go[ IFelj and the monograph jX| for further references) . 

Bose-Einstein condensates are trapped systems: A magneto-optical confinement is imposed on the 
atoms. When rotating such a system, the strength of the confinement can lead to two different behaviors. 
If the trapping potential increases quadratically with the distance from the rotation axis ('harmonic' 
trap), there exists a limiting angular velocity that one can impose to the gas. Any larger velocity would 
result in a centrifugal force stronger than the trapping force. The atoms would then be driven out 
of the trap. By contrast, a stronger confinement ('anharmonic' trap) allows in principle an arbitrary 
angular velocity. In this paper we focus on the two-dimensional GP theory for a BEG with anharmonic 
confinement. 

Theoretical and numerical arguments have been proposed in the physics literature (see, e.g., jF JS[ IFBl 
IKBp in favor of the existence of three critical speeds at which important phase transitions are expected 
to happen: 

• If the velocity D, is smaller than the first critical velocity fJci , then there arc no vortices in the 
condensate ('vortcx-frec state'); 

• If is between i7ci and ilc2, there is a hexagonal lattice of singly quantized vortices ('vortex-lattice 
state'); 

• When is taken larger than ^Ic^ , the centrifugal force becomes so important that it dips a hole in 
the center in the condensate. The annulus in which the mass is concentrated still supports a vortex 
lattice however ('vortex-lattice-plus-hole state'), until crosses the third threshold ftcs', 

• If is larger than fie.,, all vortices retreat in the central low density hole, resulting in a 'giant 
vortex' state. The central hole acts as a multiply quantized vortex with a large phase circulation. 
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In |CDYli ICYl ICRYl fR] we have studied these phase transitions using as model case a EEC in a 'flat' 
trap, that is a constant potential with hard walls. This is the 'most anharmonic' confinement one can 
imagine and serves as an approximation for potentials used in experiments. Mathematically, it has the 
advantage that the rescaling of spatial variables as e — > and/or O — > oo is avoided. The GP energy 
functional in the non-inertial rotating frame is defined as 

£^^m:=J dr {|V*p-2**0-L«' + £-2|^|4| 

where we have denoted the physical angular velocity by 257, L = —if AV is the angular momentum 
operator and B the unit two-dimensional disc. We have written the coupling constant as e~^. The 
subsequent analysis (as well as the papers [CDYll fCYl ICRYl IR] ) is concerned about the 'Thomas-Fermi' 
(or strongly interacting) limit where e — > 0. 

The simplest way to define the ground state of the system is to minimize the energy functional (jl.ll) 
under the mass constraint 

/ df 1*1^ = 1 
Jb 

with no further conditions. This is the approach that has been considered in the previous papers [CDYll 
ICYl ICRYl [R] ■ leading to Neumann boundary conditions on dB. We will refer to this situation as the 'flat 
Neumann problem' in the sequel. 

There are, however, both physical and mathematical reasons for considering also the corresponding 
problem with a Dirichlct boundary condition, i.e., requiring the wave function to vanish on the boundary 
of the unit disc. Physically, this corresponds to a hard, repelling wall which is usually a closer approxi- 
mation to real experimental situations than a 'sticky' wall modeled by a Neumann boundary condition. 
The Dirichlet boundary condition can be formally implemented by replacing the flat trap with a smooth 
confining potential of the form r'^ and taking s — >■ oo. 

Mathematically, the new boundary condition is responsible for some new aspects requiring several 
modifications of the proofs. For one thing, the density profile is no longer a monotonously increasing 
function of the radial variable and the position of the density maximum has to be precisely estimated. 
Furthermore, energy estimates have to be refined to take the boundary effect into account, and a boundary 
estimate for the GP minimizer, that was an important ingredient in the proof of the giant vortex transition 
in |CRYj ■ has to be replaced by a different approach. 

In addition to these adaptations to the new situation the present paper contains also substantial im- 
provements of results proved previously in the Neumann case. These concern in particular the uniform 
distribution of vorticity in the bulk (Theorem I l.ip and the rotational symmetry breaking (Theorem II. 6p . 
Besides, the error term in our energy estimate in Theorem 11.41 below is much smaller than the corre- 
sponding term in jCRYl Theorem 1.2]. This last improvement is due to the new method for estimating 
a potential function that we use to avoid the boundary estimate. 

From now on the minimization of (|l.ip is considered on the domain 

i^GP:={vI/ei/i(6): 11*112 = 1}, (1.2) 

where Hq{B) is the Sobolev space of complex valued functions '9 on B with /gd^'P + |V4'p) < oo and 
^'(t^ = on dB. The ground state energy is thus defined as 

E^^ := inf (1.3) 



^This limit has to be taken with care, however, because it can not be interchanged with the asymptotic limit eQ oo 
we shall consider. This point will be discussed further in ICPRYj . 
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and any corresponding minimizer is denoted by ^'^^ . This case will be referred to as the 'fiat Dirichlet 
problem'. In the following we will often use a different form of the GP functional which can be obtained 
by introducing a vector potential, i.e., 



= J d?^ I (v - tl) * ^ 

where 

A:=nAf=nrei). (1.5) 
Here (r, i?) arc two-dimensional polar coordinates and a unit vector in the angular direction. 

The GP minimizer is in general not unique because vortices can break the rotational symmetry (see 
Section [1.3p but any minimizer satisfies in the open ball the variational equation (GP equation) 

- Avl/GP - 20 • i VI/GP + 2e-2 |VI,GP|2 vl/GP = ^GP^GP^ g) 

with additional Dirichlet conditions at the boundary, i.e., 

*'^P(f) = on dB. (1.7) 
The chemical potential in (jl.6p is given by the normalization condition on ^E**^^, i.e., 

^GP^^^GP^ 1 /•d^^|^GP|4^ (1.8) 



For such a model, variational arguments have been provided in [FBj to support the following conjec- 
tures about the three critical speeds: 

f^ci (X |loge|, (1.9) 
n,, cx e-\ (1.10) 
1^03 oc e~^\\oge\-\ (1.11) 

As for the behavior of the condensate close to 0^, the centrifugal force is not strong enough for the 
specificity of the anharmonic confinement to be of importance. A consequence is that the analysis 
developed in |IM11 IIM2] (see also |AJRj for recent developments) for harmonic traps applies and leads to 
the rigorous estimate 

f^ci = |log£|(l + o(l)) (1.12) 
when e — >■ 0. In this paper we aim at providing estimates of flc2 and thus will assume that 

^ > |loge|, 

i.e., we consider angular velocities strictly above fici- The situation is then very different from that in a 
harmonic trap because of the onset of strong centrifugal forces when CI approaches ■ 

Our main results can be summarized as follows. We show that if i7 < 2{^/^^e)^^, the condensate is 
disc-shaped, while for fl > 2{y/Tr£)~^ the matter density is confined in an annulus along the boundary of 
B. In addition we prove that if 

|loge| < 0< — r, 

e-^l loge| 

there is a uniform distribution of vorticity in the bulk of the condensate. Although our estimates are not 
precise enough to show that there is a hexagonal lattice of vortices, these results support the qualitative 
picture provided in [FBj . We deduce that when e — > 

r!e, = ^(l-fo(l)). (1.13) 
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"c3< o-.,n.„,, (l + o(l))- (1.14) 



We refer to Section [TTT] for the detailed statements of these results. 

In Section fTT^ we present our results about the third critical speed. We show that if = JIq^^^I logej^^ 
with fio > 2(37r)~^, then there are no vortices in the bulk of the condensate. This provides an upper 
bound on the third critical speed 

2 

Stte^I logs I 

It should be noted right away that we do believe that this upper bound is optimal. This has been proved 
in [n] in the flat Neumann case and the adaptation of the adequate tools to the flat Dirichlct case is 
possible but beyond the scope of this paper. We hope to come back to the regime f2 oc in the future. 

We also remark that the estimates we obtain for the three critical speeds in the limit e ^ are the 
same in the flat Neumann and Dirichlct settings. In the cases of the first and second critical velocities 
this is plausible because the features that mark the onset of the transition (the first vortices and the 
appearance of the 'hole' respectively) occur far from the boundary of the trap. The independence of the 
third critical velocity of boundary conditions is less obvious but the main reason is that the maximum of 
the density is to leading order the same for both boundary conditions. 

In the regime Q > f2c3 a very natural question occurs about the distribution of vorticity in the central 
hole of low matter density: Is the phase of the condensate created by a single multiply quantized vortex 
at the center of the trap? We show that this is not the case in Section 11.31 and, as a consequence, the 
rotational symmetry is always broken at the level of the ground state, even when Q > fie,. 

Before stating our results more precisely, we want to make a comparison with the 2D Ginzburg-Landau 
(GL) theory for superconductors in applied magnetic fields (see |BBHi IFH[ ISS2| for a mathematical 
presentation). The analogies between GP and GL theories have often been pointed out in the literature, 
with the external magnetic field playing in GL theory the role of the angular velocity in GP theory. We 
stress that our results in fact enlighten significant differences between the two theories. Whereas the 
first critical speed in GP theory can be seen as the equivalent of the first critical field in GL theory, 
the second and third critical speeds have little to do with the second and third critical fields of the GL 
theory. The difference can be seen both in the order of magnitudes of these quantities as functions of e 
(which for a superconductor is the inverse of the GL parameter) and in the qualitative properties of the 
states appearing in the theories. In GP theory there is no equivalent of the normal state and there is no 
vortex-lattice-plus-hole state in GL theory. The giant vortex state of GP theory could be compared to 
the surface superconductivity state in GL theory, but the physics governing the onset of these two phases 
is quite different. The main reason for this different behavior is the combined infiuence of the centrifugal 
force and mass constraint in GP theory, two features that have no equivalent in GL theory. 

We will now state our results rigorously. The core analysis that we present below is an adaptation of 
the techniques developed in |CDY11 ICYl ICRYj for the Neumann case, but the Dirichlct condition leads 
to important novel aspects that we discuss in the sequel. 

1.1 The Regime | loge\ <^ Q <ti loge|~^: Uniform Distribution of Vorticity 

Before stating our results we need to introduce some notation. We define the density functional 

'[/]:= J^df{\Vf\'-n\^f+s-^f}, (1.15) 

for any real function /. The minimization is given by 

^GP:= inf ^GP:={/ei/i(6): / = /M|/||2 = l} (1.16) 



rGPr 
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and g is the associated minimizer (see Proposition I2.ip . In order to give a precise meaning to the 
expression 'bulk of the condensate', we introduce the following Thomas- Fermi functional, obtained by 
dropping the first term in (|1.4p or (jl.lSp : 

/ dr{p'-sWp}, (1.17) 

which is expected to provide the energy associated with the non-uniform density of the condensate. We 
refer to the Appendix for the properties of its ground state energy E"^^ and associated minimizer p^^ . 
Let us define 

y^TF ,^ supp (pTF) (1^8) 

If 51 < 2{y/Tre)^'^, A^^ ~ B. while if 51 > 2(-y7re)^^, A^^ is an annulus of outer radius 1 and inner radius 
_Rh with 1 — i?h fx (eil)"^. As we shall see below, IvJ/C^Pp is close to p^^ and thus, if 51 ^ e~^, the mass 
of ^''^^ is concentrated close to the boundary of B. 

Our result about the uniform distribution of vorticity in fact holds in a slightly smaller region than , 
namely the annulus 

Auik :={?■€ 6 : R<r< (1.19) 

where, for a certain quantity 7 := 7(e, 51) > such that 7 = o(l) as e (see Section [?751 Equation 
p.35p for its precise definition), 

^ ^ fo, if 51 < ne-\ with 51 < 

\RY,+-ie-^n-\ if 2(V^£)-i < 51<e-2|iogg|-i^ ^- > 

and i?ni is the position of the unique maximum of the density g (see Proposition [221) • It should be noted 
that R is close to i?h and i?ni is close to 1 in such a way that 

|^^^\Auik|«0(£-^f^-^) = |^^n' 

i.e., the domain ^buik tends to the support of the TF density as e — 0. Also, thanks to the above 
estimate, we have 

/ df|*GP|2 = 1-0(1), (1.21) 
i.e., the mass is concentrated in .4buik- We refer to (|2.22p . (|2.23p and (|2.32p below for precise estimates 

of i?ni- 

We now state our result about the uniform distribution of vorticity. It is the analogue of jCYl Theorem 
3.3] but here we prove that the distribution of vorticity is uniform in the whole regime | loge| ^ 51 ^ 
e~^| loge|~^ whereas in |CY| this was proved only for 51 < e^^ . 

Theorem 1.1 (Uniform distribution of vorticity). 

Let ^'^^ he any GP minimizer and e > sufficiently small. If \ loge| <C 51 ^ e~^| logej""'^, there exists a 
finite family of disjoint ball^ {Bi\ '■= {B{ri, g^)} C -4buik such that 

1. Q, < 0(51-i/2), J2 Q^ < Oin^/^) and E 9^ « (1 + £^)^\ 

2. I^GPI^ > C7(H-£51) on dBi for some C > 0. 



^Throughout the whole paper the notation B{f, q) stands for a ball of radius q centered at r, whereas B(R) is a ball with 
radius R centered at the origin. 
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Moreover, denoting by di_^ the winding number of l^**^^! ^vJ/GP qq^ introducing the measure 



Remark 1.1 (Distribution of vorticity) 

The result proven in the above Theorem implies that the vorticity measure converges after a suitable 
rescaling to the Lebesgue measure, i.e., the vorticity is uniformly distributed. However such a statement 
is meaningful only for angular velocities at most of order e~^, when the TF support A^^ can be bounded 
independently of e. On the opposite if J7 ^ e~^, A^^ shrinks and its Lebesgue measure converges to 
as e~^^l~^. To obtain an interesting statement one has therefore to allow the domain S to depend on e 
with |5| ^ as e ^ 0. 

Remark 1.2 (Conditions on S) 

We remark that the lower bound on the measure of the set 5, i.e., \S\ ^ fl""'^! log(e^rj| loge|)p, is 
important, even though not optimal, as it will be clear in the proof: In order to localize the energy 
bounds to suitable lattice cells, one has to reject a certain number of 'bad cells' where nothing can be 
said about the vorticity of VE''^^. However since the number of bad cells is much smaller than the total 
number of cells, this has no effect on the final statement provided the measure of S is much larger than 
the area of a single ceU, i.e., 0(r2^^| log(£^il| loge|)p). A similar effect occurs in |CY1 Theorem 3.3], 
where the stronger condition \S\ > C is assumed. 

Remark 1.3 (Vortex balls) 

The balls contained in the family {Bi\ are not necessarily vortex cores in the sense that each one might 
contain a large number of vortices. However the conditions stated at point 1 of the above Theorem ll.il 
have important consequences on the properties of the family. For instance, if 17 3> the last one, i.e., 
^ £)f ^ e~^ri~^, guarantees that the area covered by balls is smaller than the area of the annulus ^buik 
where the bulk of the condensate is contained. At the same time the other two conditions imply that the 
radius of any ball in the family is at most 0(17^^/^) and their number can not be too large: Assuming 
that for each ball Qi ~ f2~^/^, the second condition would yield a number of balls of order at most fi, 
which is expected to be close to the total winding number of any GP minimizer. 

An important difference between the flat Neumann and the flat Dirichlet problems can be seen directly 
from the energy asymptotics. Indeed, in the flat Neumann case (see |CY[ Theorem 3.2]) the energy is 
composed of the contribution of the TF profile (leading order) and the contribution of a regular vortex 
lattice (sublcading order). In the flat Dirichlet case the radial kinetic energy arising from the vanishing 
of the GP minimizer on dB might be larger (see Remark 1.4 below) than the contribution of the vortex 
lattice. As a result the functional (jl.l6p that includes this radial kinetic energy plays a key role in the 
energy asymptotics of the problem: 

Theorem 1.2 (Ground state energy asymptotics). 




(1.22) 



then, for any family of sets S C -4buik such that \S\^ ^\ log(e^f2| loge|)p as e — >■ 0, 




(1.23) 



As e ^ 0, 



^GP ^ ^GP ^ iog(e2r!)|(l + o(l)), 



(1.24) 



if I loge| < 51 < e^S and 



iGP 



i^°P+17|l0g£|(l + 0(l)), 



(1.25) 



ife 
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Remark I.4 (Composition of the energy) 

The leading order term in the GP energy asymptotics is given by the energy E'^^ which contains the 
kinetic contribution of the density profile (see (|1.15p ). i.e., one can decompose E'^^ as E"^^ + 0{e^^) + 
©(e^/^fJ"^/^), where the first remainder is the most relevant in the regime H, < e^^ and the second becomes 
dominant for angular velocities much larger than e~^. 

The kinetic energy of the density profile can in turn be decomposed into the energy associated with 
Dirichlet conditions cx + e^l"^^!"^ and the one due to the inhomogeneity of the profile ^ \fp^ , 
which is C(l) + ©(e^Sl^l loge|) (see Remark 2.1). The first contribution dominates for any angular 
velocity f2 <C e"^| loge|~^ and this is why it is the only one appearing in p.24p and (jl.25p . 
Note also that the kinetic energy due to Dirichlet boundary conditions is, in general, much larger than 
the vortex energy contribution, i.e., the second term in (11.241) and (jl.25p . except in the narrow regime 

e-^\\og{e''n)\-^ < O < e-^\\ogel 

where the latter becomes predominant. 

An important consequence of the above energy asymptotics is that we always have (sec Proposition 

ED 

l|l*°'l'-P^'IL^(«)-«(l)«l|/'^"b(B) (1-26) 

which allows to deduce 

/ dr IfCiPp = 1 -0(1). (1.27) 

This implies that if £7 > 2{^e)^^ ^ the mass of vp'^^ is concentrated in an annulus, marking the transition 
to the vortex-lattice-plus-hole state. We thus have 

fie. = ^(l + o(l)). (1.28) 

Note that we actually prove stronger results than (|1.2ip and (|1.27p . If fi > ^1^2 1 any GP minimizer is 
in fact exponentially small in the central hole, minus possibly a very thin layer close to r = (see 
Proposition 13. 2p . 

1.2 The Regime Q ~ e~^| log£|^^: Emergence of the Giant Vortex 

When the angular velocity reaches the asymptotic regime f2 ^ e~^| logej"-'^ a transition in the GP ground 
state takes place above a certain threshold: Vortices are expelled from the essential support of any GP 
minimizer 5''^^. The density is concentrated in a shrinking annulus where such a wave function is vortex 
free. Anticipating this transition we shall throughout this section assume that 



e^l loge| 



for some constant fip > 0. 

The bulk of the condensate has to be defined differently in this regime: We set 

Auik:={/^e6: i?> <r < l-£3/2|ioge|2| (1.30) 

where 

R> := Rh + e\\oge\-\ (1.31) 
The main result in this regime is contained in the following 
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Theorem 1.3 (Absence of vortices in the bulk). 

// the angular velocity is given by (jl.29p with Oo > 2(3tt)^^, then no GP minimizer has a zero inside 

^buik if £ is small enough. 

More precisely, for any r £ .4buik, 

<0(e-3/4|iog£|2)^^TF(^)^ (^32) 



VI;GP(^-)|2_^TF(^) 



Remark 1.5 (Bulk of the condensate) 

As the notation indicates, the domain ^buik contains the bulk of the condensate: Using the exphcit 
expression (jA.ll) of p'^^(r), one can easily verify that 

ll^^"IL^(^.„..) = l-^(|l°g^H-^), (1.33) 

which implies by (|1.32p that the same estimate holds true also for 
A consequence of this result is the estimate 

^e-,<T^, r(l + o(l)). (1.34) 

37re^| log£| 

As already noted, we believe that this upper bound is optimal, i.e., we actually have 

2 

37re2| logej 

The proof of this conjecture could use the tools of [R] but we leave this aside for the present. 



The theorem above is based on a comparison of a minimizer with a giant vortex wave function of the 
form 

/(r)exp{i(P-^)^}, 

where [ • ] stands for the integer part and G Z is some additional phase. Therefore we introduce a 
density functional 

£U[.f] [/(r)exp{z([n]-c.)^}] = 

j dr {|V/|' + ([O] - ufr-^f^ - 2(P - u:)nf^ + e-^/^} = 

df { |V/|^ + Blf^ - a\^f^ + e-^f^] , (1.35) 

where / € f^^^ is real-valued and 

B^{r):^{nr-{[n]-oj)r~^)e^. (1.36) 

We also set 

By simply testing the GP functional on a trial function of the form above, one immediately obtains the 
upper bound 

In the following Theorem we prove that the r.h.s. of the expression above gives precisely the leading order 
term in the asymptotic expansion of E'^'^ as e — and we state an estimate of the phase optimizing E^^ . 
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Theorem 1.4 (Ground state energy asymptotics and optimal phase). 

For any flo > (Stt)^^ and e small enough 

^GP ^ ^gv ^ I ioge|)-2| ioge|2). (1.39) 

Moreover E^'^ — i?5Ipt with Wopt £ N satisfying 

2 



^opt 



3>/7r£ 



(l + Odloger^). (1.40) 



Remark 1.6 (Composition of the energy) 

We refer to jCRYi Remark 1.4] for details on tlie energy E^^ (denoted E'^^ in tlrat paper). Let us 
just emphasize that in this setting the Dirichlet boundary condition is responsible for a radial kinetic 
energy contribution that was not present in the flat Neumann case and gives the leading order correction 
(X e~^/^| loge|~'^/^ to E"^^ in the asymptotic expansion of E^^ . 

A consequence of Theorem 1 1.31 is that the degree of vI/^p jg -^qW defined on any circle dB{r) of radius 
r centered at the origin, as long as 

R><r< I- e^/^\\oge\^. 

We are able to estimate this degree, proving that it is in agreement with that of the optimal giant vortex 
trial function (|1.40p : 

Theorem 1.5 (Degree of a GP minimizer). 

// r^o > 2(37r)~^ and e is small enough, 

deg{*GP,ae(r)} =f2--^(l + 0(|log£|-4), (1.41) 

for any i?> < r < 1 — £'^/^| logep. 

We note that because of the Dirichlet condition there is a small region close to dB where the density 
goes to zero. We have basically no information on the GP minimizer in this layer that could a priori 
contain vortices. The existence of this layer is the main difference between the flat Dirichlet case and 
the flat Neumann case considered in |CRY| . In particular the lack of a priori estimates on the phase 
circulation of vI/GP on dB requires new ideas in the proof. 



1.3 Rotational Symmetry Breaking 

As anticipated above, a very natural question arising from the results in Section 1.2 is that of the 
repartition of vortices in the central hole of low matter density. In particular, does one have 

*''^-.gopt(r)e'("-"-')^ 

modulo a constant phase factor, which would imply that all the vorticity is contained in a central multiply 
quantized vortex? 

We show below that this can not be the case: the GP functional is rotationally symmetric but if the 
angular velocity exceeds a certain threshold this symmetry is broken at the level of the ground state. No 
minimizer of the GP energy functional is an eigenfunction of the angular momentum, i.e. a function of 
the form /(r)e*"'' with / real and n an integer. A straightforward consequence is that there is not a 
unique minimizer but for any given minimizing function one can obtain infinitely many others by simply 
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rotating the function by an arbitrary angle. In other words as soon as the rotational symmetry is broken, 
the ground state is degenerate and its degeneracy is infinite. 

In [CDYli Proposition 2.2] we have proven that the symmetry breaking phenomenon occurs in the 
case of a bounded trap B with Neumann boundary conditions when c| loge| < ■, for some given 

constant c. We arc now going to show that such a result admits an extension to angular velocities much 
larger than £~^, i.e., the rotational symmetry is still broken even for very large angular velocities. Such 
an extension is far from obvious in view of the main result about the emergence of a giant vortex state 
discussed above: Since vortices are expelled from the essential support of the GP minimizcr, there might 
a priori be a restoration of the rotational symmetry but the behavior of any GP minimizcr inside the 
hole B{Rh) remains unknown. 

Theorem 1.6 (Rotational symmetry breaking). 

// e is small enough and eCt large enough, no minimizer of the GP energy functional (jl.ip is an eigen- 
function of the angular momentum. 

We note that it is proved in jAJR] for a related model that the ground state is rotationally symmetric 
if < fici and e is small enough. Theorem 11.61 shows that the symmetry, broken due to the nucleation 
of vortices, never reappears, even when fl > flc^- 

1.4 Organization of the Paper 

The paper is organized as follows. Section [2] is devoted to general estimates that will be used throughout 
the paper. We then prove our results about the regime |loge| <^ fl <^ e~^|loge|~^ in Section [31 The 
analysis of the energy functional (jl.l5p is the main new ingredient with respect to the method of [CY| . 
We adapt the techniques developed in that paper for the evaluation of the energy of a trial function con- 
taining a regular lattice of vortices. The corresponding lower bound is proved via a localization method 
allowing to appeal to results from GL theory |SS1[ ISS2j . The inhomogeneity of the density profile is dealt 
with using a Riemann sum approximation. 

Section U is devoted to the giant vortex regime. Our main tools are the techniques of vortex ball con- 
struction and jacobian estimates, originating in the papers [Sal [31 US] (see also |SS2| '). We implement 
this approach using a cell decomposition as in |CRYj . New ideas are necessary to control the behavior of 
GP minimizers on dB. 

The symmetry breaking result is proved in SectionjSl Following jSeirj . given a candidate rotationally sym- 
metric minimizer, we explicitly construct a wave function giving a lower energy. Finally the Appendix 
gathers important but technical results about the TF functional and the third critical speed. 

2 Preliminary Estimates: The Density Profile with Dirichlet 
Boundary Conditions 

This section is devoted to the proof of estimates which will prove to be very useful in the rest of the paper 
but are independent of the main results. We mainly investigate the properties of the density profile which 
captures the main traits of the modulus of the GP minimizer |^'^^|: More precisely we study in details 
the minimization of the density functional £'^^ p.l5p and prove bounds on its ground state energy i?^^ 
(|1.16p and associated minimizers g. 

The leading order term in the ground state energy E'^^ is given by the infimum of the TF functional 

™' i?^" inf := {P ^ L\B) : p > 0, = l} . (2.1) 
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We postpone the discussion of the properties of E as well as the corresponding minimizer p to the 
Appendix. 

Proposition 2.1 (Minimization of £'^^). 
//O < loge|"2 as e ->-0, 

j^TF < ^GP < ^TF ^ 0(^-1) ^ (el/2f^3/2j (2.2) 

Moreover there exists a minimizer g that is unique up to a sign, radial and can he chosen to he positive 
away from the boundary dB. It solves inside B the variational equation 

- Ag - n^r^g + 2£- V ^ ^GP^^^ ^2.3) 
with boundary condition g{l) = and ffi^ ~ E'*^^ + Il5|l4- 
Remark 2.1 ( Composition of the energy 

The remainders appearing on the r.h.s. of p.2[) can be interpreted as the kinetic energy due to Dirichlet 
boundary conditions: The bending of the TF density close to r = 1 in order to fulfill the boundary 
condition produces some kinetic energy which is not negligible and can be estimated by means of the 
trial function used in the proof of the above proposition, i.e., 0{e~^) as long as f2 < e""'^, and 0{e^^'^fl^^^) 
for larger angular velocities. Note indeed that the second correction becomes relevant only if f2 > e~^. 
The orders of those corrections can be explained as follows: li ft < the TF density goes from its 
maximum of order 1 to in a layer of thickness ~ e (because of the nonlinear term) , yielding a gradient 
^ e^^ and thus a kinetic energy of order e~^. If $7 ^ the thickness of the annulus where g varies 
from y/eTi to becomes of order e^/^fj^^/^ and the associated kinetic energy is 0{e^^^Q'^^^). 
Note that in both cases the kinetic energy associated with the boundary conditions is much larger than 
the radial kinetic energy of the profile \/ p^^ which is 0(1) in the first case and ©(e^fi^j logej) in the 
second one (see |CY[ Section 4]): The condition il ^ e^^ \ logej"^ is precisely due to the comparison of 
such energies for large angular velocities. 

Finally wc point out that, if 17 ^ £~^, the correction of order due to Dirichlet boundary conditions can 
become much larger than two terms of order Vl^ and e^fJ'' contained inside E"^^ (see the explicit expression 
(jA.3p in the Appendix), so that the upper bound could be stated in that case E'^^ < tt^^e^^ +0{e~^). 



Proof of Proposition 

The lower bound is trivial since it is sufficient to neglect the positive kinetic energy to get E'~^^ > E'^^ . 
The upper bound is obtained by evaluating £'^'^ on a trial function of the form 

ft.Ur) = c^^^nir) (2.4) 

where c is the normalization constant and < (r) < 1 a cut-off function equal to 1 everywhere except 
in the radial layer [1 — (5, 1], 5 <C (1 + e^l)~^, where it goes smoothly to 0, so that / satisfies Dirichlet 
boundary conditions. The density p{r) coincides with p"'"^(r) if il is below the threshold 2{y/TTe)~'^ and 
is given by a regularization of p'^^ above it, i.e., if efi > 2/y^, we set as in jCYI Eq. (4.9)] 

{0, if < r < i?h, 

f72pTF(^j^ ^ _ ^j^)2^ if i?h < r < i?h + n-\ (2.5) 

p^^{r), otherwise. 

Notice that p differs from p"^^ only inside the interval [i?h, i?h + f^^^] and 



p{r)=p'^''{r)+0{e^n). 



(2.6) 
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In order to estimate the normalization constant we use the bound p^^ < C(l + efi), which imphes 

c-^ = f drpil> [ drp^^-2TT f dr r{l - £,l)p - Ce^ > 1 - C[{en + 1)S + e^]. (2.7) 
Jb Jb Jl-S 



IB JB Jl-S 

The kinetic energy of /trial is bounded as follows: 



^dr- |V/tHai|' < 2c' ^df + p\W^d\^} < C[e'n^\ loge| + {en + l)S-^], (2.8) 

where we refer to |CY1 Eqs. (4.14) and (4.15)] for the estimate of the kinetic energy of p. 
The interaction term can be easily estimated as 

i / drfti. < ' + + Vll < Wp^II + + e-rs. (2.9) 



e 

To evaluate the centrifugal term we act as in |CY[ Eqs. (4.44) - (4.46)]: With analogous notation 
-n^ f d7^rVtHai--27rl)'+47rl]2 f\rr T dr' r' /trial (^') < 

JB Jo JQ 

- 27rl72 + AttH"^ [ drr [ dr' r'p^^{r') + C[n + {e^ + ^-^y^S] = 
Jo Jo 



/ drr^p^^ + C{n + e-'^S + n^d), (2.10) 



8 

where we have integrated by parts twice and used (j2.6p . (|2.7p and the normalization of /trial- Hence one 
finally obtains 

S^'' [/trial] < + C[e^n^\ logel + n + {en + 1)5-^ + + e-^ fS]. (2.11) 

It only remains to optimize w.r.t. i5, which yields 5 — e, if il < e^^ , and S — e^^^Q^^^^ otherwise, and 
thus the result. □ 

A crucial property of the density g is stated in the following 

Proposition 2.2 (Behavior of g). 

The density g admits a unique maximum at some point < i?m < 1 • 

Proof. The method is very similar to what is used in |CRY[ Lemma 2.1], although in that case one 
considers the Neumann problem. After a variable transformation — > s the functional S'^^ becomes 

TT ds {s|V/|2 - n'sf + e-V'} , (2.12) 
Jo 

and the normalization condition ^ 

/ dsg^ = TT-^. (2.13) 

"'0 

We first observe that the Dirichlet boundary condition implies that g cannot be constant, otherwise we 
would have g = everywhere, contradicting the mass constraint. 

Suppose now that g has more than one local maximum. Then it has a local minimum at some point 
s = S2 with < S2 < 1, on the right side of a local maximum at the position s = si, i.e., si < S2- For 
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< e < g^(si) — g^{s2), we consider the set = {0 < s < S2 : 5^(si) — e < .9^(s) < g^{si)}: Since g is 
continuous, the function 

^e):^J^dsg^ (2.14) 

is strictly positive and $(e) — as e — 0. Likewise, for any k > 0, we set J7k = {si < s < 1 : 5^(52) < 
S^(s) < g^(s2) + k}, so that 

r(K) := / dsg^ (2.15) 

has the same properties as 

Hence, by the continuity of g, there always exist e, k > 0, such that g^(s2)+K < _g^(si) — e and $(e) = r(K). 
Note that this implies that and are disjoint. 
We now define a new normalized function g by 

rg2(5i)2_e, if,sel^, 
f{s):=lgHs2) + K, ifsej-, (2.16) 
[(7^(s), otherwise. 

The gradient of g vanishes in the intervals and and equals the gradient of g everywhere else, so that 
the kinetic energy of g is smaller or equal to the one of g. The centrifugal term is lowered by g, because 
— s is strictly decreasing and the value of g^ on le is larger than on J7k- Finally since mass is rearranged 
from to Jtj, where the density is lower, \\g\\'^ < WgW^- 

Therefore the functional evaluated on g is strictly smaller than E'^^ , which contradicts the assumption 
that g is a minimizer. Hence g has only one maximum. □ 

The energy asymptotics ()2.2p implies that the density g^ is close to the TF minimizer p^^: 

Proposition 2.3 (Preliminary estimates of g). 



lL2(e) 



<0(ei/2+e5/4f73/4)^ (2.17) 



9\Rn.) = Ml^^r,) < l|P^nL=^(e) {l + OiV-e) + Oie^^'n'/')) . (2.18) 
Proof. See, e.g., [CRYI Proposition 2.1]. Note that p.l7|) implies the bound 

I^GP _ ^TFI < ^ ^ ^)l/2 / ^-3/2 ^ g-3/4^3/4\ ^ (2.19) 



which yields ^J^. □ 

Next proposition is going to be crucial in the proof of the main results since it allows to replace the 
density g^ with the TF density p^^: On the one hand, using the fact that the latter is explicit, this result 
will be used to obtain a suitable lower bound on g'^ in some region far from the boundary and, on the 
other hand, it implies that the boundary layer where g goes to is very small. 

Proposition 2.4 (Pointwise estimate of g). 

Ifn< Cle-'^ with n < 2/ VtF ase^O, 

\g^r)^p'^^{r)\<0{V~e), (2.20) 
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for any < r < 1 — 0{e\ loge|). 

On the other han^ if 2{y/l^e)-^ <^<£^^, 

\g^{r) - pTF(r)| < ©(e^/^f^s/^), (2.21) 
for any i?h + 0{e-^n-^\ log£|-2) < r < 1 - 0{e^/^n-^/^\ log£|3/2). 
Remark 2.2 (Position of the maximum of g) 

The pointwise estimates (j2.20p and (|2.2ip give some information about the position of the maximum 
point of g. Assuming that < e^^, one has the lower bound 

ff'(i?m)>P^^(l-£|loge|)-Cv/i, 

since (j2.20p holds true up to a distance e\ logej from the boundary. Hence one immediately obtains 

R„,>l-Oie-^^^n-^), (2.22) 

which for e~^^^ ^ Q < s^^ implies that i?„i 1 — o(l). For smaller angular velocities the above inequal- 
ity becomes useless: Since p^^ is approximately constant in those regimes, i.e., p'^^(r) — tt^^ + 0{e^n'^), 
the pointwise estimate (j2.20p is too rough to extract information about the maximum of g. 
On the opposite if 2{^/^^e)^^ ^ ^ ^ we get from (|2.2ip the following: Either i?ni > 1 — 

e-^/^fi"-'^/^! logep/2 or the pointwise estimate (|2.2ip applies at i?m yielding g^{Rm_) < p^^iRm) + 
CeV^n^/'i and 

<?'(i?m) > P^^(l - ei/2o-i/2| log£|3/2) - Cey^n'^\ 

so that, in any case, 

i?m > 1 - Oie-^/^n-^/^), (2.23) 
since ei/2f^-i/2| iog^|3/2 ^ e-^/^Q-^/^ 

Remark 2.3 (Improved pointwise estimates of g) 

Thanks to the remark above, it is possible to refine the estimates (|2.20p and (|2.2ip and extend them up 
to the maximum point of g. More precisely one has the following: If < Cle^^ with f2 < 2/^7?, 

\g{r)~p^^{r)\<0{^e), (2.24) 

for any < r < max[i?„i, 1 - e\ loge|]. If 2(0Fe)-i <n< e^^ 

\g{r) - pTF(^)| < o(e7/4^5/4)^ (2.25) 

for any R^ + e-^Vt~^\\oge\-'^ <r< max[i?m, 1 - e^/^n-^/^\\oge\^/^]. 

The extension can be easily done in the first case {il < Cte~^) by noticing that one can suppose R„^ > 
1 — e| logej (otherwise the bound is given by the original result), so that (|2.24p follows from (|2.18p together 
with 

and the fact that g is increasing in i3(i?ni)- 

In the other regime the key point is the estimate (|2.23p . which implies 

max |pTF(l) _^TF(^) I < ^^7/4^5/4^ 

l-e-i/*0-3/-i<r<l 



^This second estimate applies also if Q = 2{y/ne) ^{1 — o(l)), in which case has to be set equal to 0. 
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Proof of Proposition [ 

The proof is done exactly as in jCRYl Proposition 2.6], so we highlight only the main differences. 
The result is obtained by exhibiting suitable local sub- and super-solutions to the variational equation 



where the function p is given by 

By (EUl), iin< 2(v^e)-\ 
whereas if > 2(v^e)~^, 



-Ag = 2s-^p-g^)g, (2.26) 
pWi^H^'A^'^+^W). (2.27) 



IL^(e)^^Vi, (2.28) 

|/5(r)-pTF(^)| <C'e7/4^5/4^ (2.29) 
for any r > _Rh- 

In order to apply the maximum principle one needs a lower bound on the function p in the domain 
under consideration and it is provided by the above estimates: In the fist case, i.e., if $7 < fie^^, 
p'^^ir) > C{Cl) > and the pointwise estimate (|2.28|) guarantees the positivity of p everywhere; otherwise, 
if 2(0Fe)-i <n< £-2^ the condition r > i?h + 0{e-^n-^\\oge\-^) yields 

pTF(r) > Cen\\oge\-^ > 0{e^/^n''^^) > \p{r) - p^^ir)\ , 

so that p{r) > Cen\ loge|^^ > in the region considered. 

The rest of the proof is done as in jCRYl Proposition 2.6] in a local annulus [ro — S, ro + S] with S = e\ loge|, 
if f2 < Qe~^ , and S = e^^'^fl~^^'^\ logep''^ otherwise. Note that the lack of monotonicity of the density 
profile g prevents a straightforward extension of the estimate to the whole support of p^^ . □ 

For angular velocities larger than the threshold 2{y/Tre)~^ the TP density develops a hole centered at 
the origin of radius i?h (see the Appendix) and in this case one can show that the density g is exponentially 
small there: 

Proposition 2.5 (Exponential smallness of g inside the hole). 

Ifn> 2(0Fe)-i ase~^0, 



Hi 



g'{r) < Can exp ^ - j— ^ } (2-30) 



for any f G B. Moreover, if ^ > 2{y/TTe) ^ + 0(1), there exist a strictly positive constant c such that for 
any f such that r < i?h — 0(e^/^), 

g^r)<Csne^p[--^}. (2.31) 

Proof. See |CRY[ Proposition 2.2]. Note that in the proof of the second estimate, the condition il > 
2(0re)"^ + 0(1) is needed in order to guarantee that i?h > 0{s''^'^). □ 

The pointwise estimates (|2.24p and (|2.25p and the exponential smallness stated in the proposition 
above have some important consequences as, e.g., an improved estimate on the density g close to the 
boundary of the trap: 

Proposition 2.6 (Estimate of i?m and estimate of g). 

Ife-^ < f7 < s-^, 

Rn.>l- Oie-'/'n-y^), Il5lli.(f,\6(«„,)) ^ Oie'^'n'/'). (2.32) 
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Proof. Without loss of generality we can assume ft > 2{y^e) ^, since the proof in the other case, i.e., 
without the hole, is even simpler. Because of the normalization of both p^^ and g, we have 

i3\6(fl„) JB{fl.,„) 

The monotonicity of g in B{R„,) and the bound fTI^ yield g'^iRh) < 5^(i?h + e~^l^"^| loge|-i) < 
ef2|loge|~^, so that, setting for convenience Rq := Rh + e~^^~^\loge\~-^ and using the exponential 
smallness p.3ip . one obtains 

/ df|pTF_92| < /■ drg' + Ce'/'^n~'\\oger'+ f dr\p^^-g'\< 

JB(Ro) JB{Ri,-eVa) Je(i?o)\e(/?h) 

C ||pTF _ ^2||^^^^^ e-i/2f^-i/2| iog£|-i/2 + ioge|-i < Ce3/4f]i/4| ioge|-i/2. 

For r > Rq one can apply the pointwise estimate (j2.25p . which yields 



JB(R,^)\B(Rn) 



I B{R„,)\B{Ro) 

Collecting the above estimates one therefore has 

/ df(pTF_52) < ^^3/4^1/4^ (2.33) 

JB\B{Rr„) 

On the other hand by (l^?^ . g^{Rm) < p^^{Rni) + Ce'^/^n^/^), so that 

/ dr{p^^ ^g^)>'^ f df (r^ - Ri) Cey'n^'\l - O > 

JB\B{R^) ^ JB\B{Rrr,) 

\ne'n'{l-Rir-Cey^n'/\l-R^,,), 

which gives the estimate of R^. 

Since the argument leading to (j2.33p is symmetric in g"^ and p"^^ , it is also true that 

/ drg'< f dr + C£3/4r!V4 < ^^3/8^1/8^ (2.34) 

JB\B{Rra) JBXBiRrn) 

due to the lower bound on i?i„ (|2.32l) . □ 

3 The Regime | log£| <C <C loge|~-^ 

This section contains the proof of the main results stated in the Introduction for the regime |loge| ^ 
f2 <C e~^| loge|~^. We also prove some additional estimates, which are basically corollaries of the main 
results and will be used also in the analysis of the giant vortex regime. 

3.1 GP Energy Asymptotics 

The most important result proven in this section is the GP ground state energy asymptotics: 
Proof of Theorem 

The result is proven by exhibiting upper and lower bounds for the GP ground state energy. 
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Step 1. For the upper bound we evaluate the GP functional on the trial function 



(3.1) 



where c is a normalization constant. $(r) the phase factor introduced in ^CY, Eq. (4.6)] and ^ a cut-off 
function: More precisely, using the complex notation C, ^ x -\- iy E C for points f = {x,y) G M^, we can 
express $ as 

c-c. 



w) ■■= n 



ic-c.r 



(3.2) 



where we denote by £ a finite, regular lattice (triangular, rectangular or hexagonal) of points fi G B such 
that the corresponding cell is contained in B: Each lattice point r,; lies at the center of a lattice cell 
and the lattice constant i is chosen so that the area of the fundamental cell Q is 



(3.3) 



Thus £ = Cfi and the total number of lattice points in the unit disc is $7(1 — 0{Vl ^^'^))- In order to 
get rid of the singularities of the phase factor $ at lattice points, we define the function 



1, if |c-Cd>i, vce/:, 

t-^\C - Q\, if IC - Q\ < t, for some Q&C 



(3.4) 



where t is a variational parameter satisfying the conditions min 

The normalization constant takes into account the effect of the cut-off function ^ and it is not difficult 

to see that 1 < < 1 + Cnt'^. 

We start by computing the kinetic energy of ^'triai: 



dr 



V-iA] Atrial 



dr |V(50r 



(3.5) 



The first term in the expression above can be estimated as follows: 

/ dr |V (5^)1^ - / df |V5|^ < ic^ / dfVg^ . V^^ + C (f7 + Qe-h^ + e^/^il'^/H^) < 

JB JB JB ^ ' 

P,^JdB{fi,t) JB ^ ' 



where we have used the bounds < t-\ \A^^\ < Ct-'^ and ||Vg||2 < C{£-^ + £1/2^73/2) (g^^, 
We have also used the fact that 

V sup g\r)<Cn, 

which can be seen as a consequence of the upper bound < C{£^ + 1) in addition to the exponen- 
tial smallness (|2.3ip . which allows to estimate the above quantity as the number of cells contained in 
supp {p^^) times eVL -f 1, i.e., 0(r2). 

In order to estimate the last term in p.5[) . we act exactly as in |CY[ Proposition 4.1]. The estimate (4.37) 
in |CYj ■ that is obtained by making use of an analogy with an electrostatic problem, reads in our case 



2„2 



< (1 + Cm^'-') sup g\r) {ttI \og{t'n)\ + 0{1)) 



(3.7) 



fiec 



Critical Speeds in the GP Tlieory - CPRY - Mardi 10th, 2011 



19 



It remains then to use the Riemann sum approximation and the normahzation of to estimate the sum 
in the above expression: If $7 < Cle~^ for some < ^/y/ir, we can simply use (|2.24[) to replace with 
p^^ and proceed as in the proof of Proposition 4.1 in |CY| . obtaining 



df^2g2 



< 



1 + o{m'/^) + oiV^)] \Q\'^ (ttI \og{t^n)\ + o(i)) 



(3.8) 



Note that inside each cell supp'^^ — inf p^^ < Ce^fl^^^ <^ ^/s, so this error term can be absorbed in the 

0{y/e) in the equation above. 

In the opposite case, if i7 > 2(v^£)~^, we set 



with 

so that ^ - i?h < e~^fl^^ and 



V := {fe B: r>R}, 
R:=R^ + e-^n-^\\og{e^n\ \oge\)\-\ 



P^^ir) > ^en\loge\-\l-oil)), WeV, 



(3.9) 
(3.10) 

(3.11) 



since | log{e'^n\ loge|)| < | loge|(l + o(l)). 

Now we can replace g^ with p^^ inside V by means of (|2.25|) . Moreover in the region r < Rwe can use the 
exponential smallness (|2.3ip . if r < i?h — e^^^, and the pointwise bound g^{r) < g'^{R) < Cefl\ logel"""^, 
if Rh — e^^^ < r < R, which follows from (|2.25p and the monotonicity of g'^{r) in B{Rm). The result is 
the upper estimate 



2„2 



< 



l + C(tOi/2+e3/4^i/4|iQgg|)j ^ sn^p^^{r){n\\og{eVL)\+0{l)) 



CeVl\\oge\-^ {r: Rh - e''^^ < r < R} [Qf^ | loge] + Cel^l loge] exp j-ce-^/^} < 
l + C(tr!i/2+e3/4^i/4|jQg^|)l SUV p^^{r){7T\log{t^n)\+0{l)) +Cn\\oge\-^ < 



1 + C{m^'^ + e3/4i]i/4| iog£| + £^^1/2)] iog(t2f7)| + o{l)) + CQ.\ logepi, (3.12) 

where we have used the estimate supp"^^ — inf p"""^ < Ce^Vt^^^ inside any cell C V. 

The estimate of £:'^^[|*triaifl can be obtained as in [CTl Eqs. (4.42) and (4.48)] (see also ((^1^ ): 



£ 



TF 



1^- 



trial 



< (1 +Cm2)£ 



drg'^ - I drr^g^ + Ce~'nH''' < 



-lo2+2 



f [^2] + c [e-^nt^ + e-'nh^] . (3.13) 



To conclude the proof of the upper bound it only remains to choose the variational parameter t: In 
the regime ft < 0.e~^, f2 < i/y/ir, we take t = e so that the remainder occurring in the above estimate 
becomes 0(51) as in (|3.6p and (|3.8p . whereas, if f2 > 2(-y/7re)~^, the remainder in (j3.13p leads to P = eil.~^ 
in order to recover the same error term 0{fl) as in (|3.6|) . In p.Sp there is an additional remainder of 
order 0{efl^^^\ loge|) which might become larger than f2 for very large angular velocities and is due to 
the Riemann sum approximation. 
The final result is therefore 



f [Atrial] < E"^"- + n\ iog{s'n)\ + OiQ), 



(3.14) 
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if 1 < f2 < e"\ and 



f [Atrial] < + n\ logel + 0{n) + log£|), 



(3.15) 



if 0<e-2. 

Step 2. The starting point of the lower bound proof is a decoupling of the energy which can be 
obtained by defining a function u{r) as 



u{r) := g-\r)^'^^{r). 



(3.16) 



Note that, thanks to the positivity of g, the function u is well defined in the open ball {r : r < 1}. 
By means of this definition and the variational equation (|2.3p , one can use a method originating in |LM) 
to decouple the energy (see, e.g., }CRY[ Proposition 3.1] or [Sel Lemma 2.2] where a Dirichlet boundary 
condition also appears) to obtain, using the normalization of both vj/'^P and g. 



We deduce the lower bound 

^GP^^GP[^GP]>^GP 

by restricting the last integral to 2?, with 
V := 



drg 



B 



(V - iA)u 



+ £-V (1-1^1- 



Afg 



{V-iA)u +e-2g2 |y|2) 



~ _ f{f e B: r< 1 -e| logel}, if f7 < Vle-^ , with n < 2/V^, 
\BiRn,)r\V if f7 > 2(0Fe)-i. 



The pointwise estimates (|2.24p and (|2.25p allow the replacement of g^ with p^^ : 

2 



> ^GP 



1-C V^ + £3/4^71/4 



dfp 



TP 



V 



(V - iA)u 



e-V^(l-|un^ 



Moreover as in |CY1 Section 5] we define another regular (square) lattice 

C := 1?^, = (m£, ni), m,ne1: Q' C p| , 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



where is the cell centered at and the lattice spacing satisfies the same conditions as in jCYl Eq. 
(5.16)], i.e., 

I log£|i/2r!-i/2 < I < min [1, {en)-^\\og{e^Vl\ loge|)|-i] , (3.21) 

so that 



sup y {r) ' p'' {r,)\ <CeVt£\\oge\p''{n). 



Hence p.l9p yields the lower bound 

\ - 0{^) - 0{e^/^n^^^\\oge\) - 0{en£\\oge\)\ ^ p^^{r,)£'^'^[u] > 

(l-o(l))^pT^(r,)f('^)M, (3.22) 



^GP _ ^GP > 
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where S'^'^ is defined as in [CYl Eq. (5.18)], i.e., 



(V - iA)u\^ + e-^p^^in) (1 - . (3.23) 



After a suitable scaling the energy above can be seen as a Ginzburg-Landau energy with a fixed external 
field /lex in the range | loge| ^ /lex ^ where e is a new small parameter. We can thus use the lower 
bound for the Ginzburg-Landau energy (see [SSli ISS2| ) as in |CY[ Proposition 5.1]. The result is 

8^'^[u] >17?|log(min[£,£2l7])|(i_o(l)), (3.24) 

for any | logej ^ f2 <C £~^| loge|^^ 

To complete the proof if suffices then to use, for any Q, ^ e^^, the estimate p^^{r) > tt^^{1 — o(l)), 
which yields 

J2 P^"^^^) ^ (1 - - log£|))|S||Qri > (1 - o{l))t\ (3.25) 

and thus the result. On the other hand, if S7 > , a simple computation (see, e.g., (|2.34p ) using the 
estimates ^(TH^ and ([O^ gives 

which implies 

E /^^^(^O > (1 - o{l))t^ [ df pTP(r) > (1 - o(l))i-2, (3.27) 



thanks to the normalization of p^^ . 

By putting together (|3.22p . (|3.24p . p.25p and p.27p . one obtains the lower bound matching (|3.14p and 
3.2 Estimates for GP Minimizers 

The GP energy asymptotics has many important consequences on the asymptotic behavior of GP min- 
imizers: For instance the upper bounds p.l4p and p.l5p immediately imply the convergence of any 
minimizing density [vJ/^Pp to the TF density p^^: 

Proposition 3.1 {L^ convegence of |^''^^p). 

As e 0, if I logel < < 

111*°^' - P^^IL^(a) < Oie'/') + 0{sn'/'\\ogie'n)\'/'), (3.28) 
whereas, if e^^ <C O <C e^'^ , 

Proof See [CRYl Proposition 2.1]. □ 

Acting as in the derivation of jCRYl Eq. (2.8)], one can show that the above estimates imply a 
bound on the chemical potential /j,*^^ occurring in the variational equation (|1.6p solved by ^'^^i 

I GP _ ^TFI < 0(^-3/2) ^ 0(e-l/2f^l/2| log(g2^)|l/2)^ (3 gg) 
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if I logel < < e"\ while, for < f2 < e"^, 

I^GP _ /,TF| < 0(^^2) ^ 0(£-l/2^| loge|l/2). (3.31) 

Such estimates can in turn be used to prove a pointwise upper bound for |v]>gp |2 (ggg [CRYl Proposition 
2.1]), i.e., 

||^GP||2 < TF.^^ fl + Oie'/^)+0{sy'n^/'\\og{s'n)\'/^), if |loge|«17<.-i, 

II ^ ^ \l + 0(£2l7)+0(£i/2|iog£|i/2), ife-i«0«e-2. ^ 

We finally state another very useful pointwise estimate of vpGP analogous to |CRYI Proposition 2.2] 
and Proposition [5?5] As is the case for the density profile g, if the angular velocity is above the threshold 
2(Y^e)~^, any GP minimizer is exponentially small inside the hole B{Rh)- 

Proposition 3.2 (Exponential smallness of 'I''^^ inside the hole). 

If ^ ^ + C'(l); as e ^ 0, there exists a strictly positive constant c such that for any f such 

that r < i?h - ©(e^/^); 

|^^P(r)|' < exp{--^} . (3.33) 
Proof See [CRYl Proposition 2.2]. □ 

3.3 Distribution of Vorticity 

We are now able to prove the uniform distribution of vorticity: 
Proof of Theorem 

The proof follows very closely the proof of jCY[ Theorem 3.3] and relies essentially on |SS1[ Proposition 
5.1]. 

The argument has to be slightly adapted depending on the value of the angular velocity: For any 
< f2e~^, with Ct < 2/y^, the proof of |CYi Theorem 3.3] applies with only one minor modification, 
since the cells in the lattice £ occurring in the lower bound proof do not cover the whole of B. However, 
since the region covered by cells tends to A^^ as e — )- and the area of the excluded set close to the 
boundary is of order O(e|loge|), i.e., much smaller than the cell area, such a difference in the lattice 
choice has no consequences for the final result. 

We now discuss the modifications in the regime ^ <C e^'^l logej""'^ which was not taken into 
account in [CYl Theorem 3.3]. The starting point is the localization of the energy bounds p.lSp . (|3.22p 
and (|3.24p . which can be rewritten as 

^ pTF(^^) -f}P|i^g^|| < ^f^pjioge] ^ p-^^in), (3.34) 

for some 

77 = 77(e, ri) < 1 as e 0. 

In order to obtain a similar estimate inside one lattice cell, one first needs a suitable lower bound on the 
density p^^ and this can be obtained by restricting the analysis to the bulk of the condensate, i.e., 

Aulk :={re6: R < r < Rn^} 

where, if » e^^, ^ is given by 

R :^ Ri, + ls-'n-\ 7:= Jlogr,]-!. (3.35) 
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We then have, for some C > 0, 

P^^(r) > Cen\ log 771-1 on Auik- (3.36) 

Moreover, the localization of the energy estimate requires that a certain number of bad cells be rejected: 
As in [CYI Theorem 3.3] we first introduce a new small parameter 



so that I \og{e^n)\ = I loge|(l + o(l)) and ((334| yields 



2£|log£| 

' ^ ' < 1, (3.37) 



5^ pTF(^^)^W[^] < (1 ^ pTF^r,)nP\logie'n)l (3.38) 



where 



£i'^[u] := / df 



(V-a)?/ + e-2(i_|u|2)2i (3 39) 



with r]{e, O) as e 0. 

We then say that a cell C .4buik is a good ceU if 

S^iu] < (1 + V^)f]P| iog(e2^)|, (3.40) 
whereas the cell is bad if the inequality is reversed. 

Now given any set S C -4buik such that \S\ ^ the upper bound (j3.38p . the definition of bad cells, 
p.36p and the upper bound p^^ < 0{en) imply that 

NB<Vvi''N^^\\og7j\N<$:N, (3.41) 

where Nb and A'' stand for the number of bad cells and the total number of cells contained inside S 
respectively. 

On the other hand by the definition (|3.40p good cells satisfy the assumptions of |SS1[ Proposition 5.1] 
and therefore one can construct a finite collection of disjoint balls {Bi} :— {B{fi, gi)} such that \u\ > 1/2 
on the boundary of each ball and Qi < 0{n^^^^). Hence one can define the winding number di^^ of u on 
dBi, which coincides with the winding number of ^1''^^ and, using |SS1[ Proposition 5.1] 

27r^d,,e 0|2(l + o(l)), 27r^ ^^^"^(1 +o(l)). (3.42) 

The rest of the statement of Theorem 11.11 easily follows by noticing that one can always take £ = 
n~i/2| log(e2j7| loge])], which satisfies p.2ip . obtaining the lower condition on the area of the set S. □ 



4 The Giant Vortex Regime Q ~ £ ^|loge| ^ 

As a preparation for the proof of the main results contained in Theorems 11.31 and 11.41 we formulate and 
prove in Section [4.11 some important propositions about the properties of the giant vortex density profiles. 
The proof of the absence of vortices in the bulk will follow the analysis of the ground state energy asymp- 
totics, which is achieved in several steps. The main ingredients are the energy decoupling (Section I4.2p . 
the vortex ball construction and the jacobian estimate (Section |4.4[) . Each individual step is analogous to 
the corresponding one contained in |CRYj and we will often omit some details, only stressing the major 
differences with the analysis of |CRYj and referring to that paper for further details. 
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4.1 Giant Vortex Density Profiles 

In this section wc investigate the properties of the giant vortex profiles and the associated energy func- 
tional defined in (jl.35p . Actually for technical reasons which will be clearer later we consider a functional 
identical to (|1.35p but on a different integration domain, i.e., 

A := {f e B: r> i?<}, (4.1) 

where i?< < i?h is suitably chosen in order to apply some estimates: All the conditions on i?< occurring 
in the subsequent proofs are satisfied if we take 

i?< Rh - e^/\ (4.2) 

More precisely we define 

^""^ := {/ e H\A) :/ = /*, = 1, / = on Ob] (4.3) 

and set, for any 

J^df{\Vf\'+Blf - . (4.4) 

We recall that 

The associated ground state energy is 

^r:= inf £^'[f] (4.5) 
and we denote guj any associated minimizer. 

The TF-like functional obtained from (j4.4p by dropping the kinetic term is denoted by (see (|A.4p ) 
and its minimization discussed in the Appendix. 

Proposition 4.1 (Minimization of £^^). 

Ifn(x £-2|iog£|-i and \uj\ < 0{e-^/^\\oge\-^/^) as e ^ 0, then 

^TF < ^TF < ^gv < ^TF ^ 0(£-5/2| loge|-3/2) < ^TF ^ 0(e-5/2 | j^g £| -3/2) ^ (4 g) 

There exists a minimizer that is unique up to a sign, radial and can be chosen to he positive away 
from the boundary dB. It solves inside A the variational equation 

- + Blg^ - n\^g^ + 2e-^gl = fi^Jg^, (4.7) 

with boundary conditions gu{l) = and gi,{R<) = and fi^ = E'g^ + e^^ ||,9w|l4- 
Moreover g^^ has a unique global maximum at with i?< < < f . 

Remark 4-1 (Composition of the energy) 

Unlike the flat Neumann case, the remainder in the r.h.s. of (j4.6p is of the same order even if the refined 
TF energy E^^ is extracted. The reason is that such a remainder is actually due to the radial kinetic 
energy of the giant vortex density profile and in particular to Dirichlet boundary conditions. 
In order to give some heuristics to explain such a difference with the flat Neumann case, it is indeed 
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sufficient to note that, by the pointwise estimate (j2.2ip . the density g^^ goes from its maximum value 
~ ei/2^|i/2 ^e-i/2|ioge|-i/2 to in a region of width at most 0{e^^^n-'^/^\\ogs\^/^) = 0{e^/^\loge\^). 
This yields an estimate for the kinetic energy of in that region as 0{e~^^'^\ logej"''), i.e., approximately 
the same remainder as in ()4.6p . which is in any case much larger than the difference between the TF 
energies E'^^ - E^^ (see (|A.5|) ). 



Proof of Proposition \4-l\ 

The proof of Proposition 12.11 applies to the functional noticing that 

\BUr)\ <nsup{r-^-r)+C\uj\<C{e-^ + \uj\), (4.8) 

which implies that the term in the functional (see the second expression in (|4.4p ) is always smaller 
than the remainder in (|2.2p . provided |a;| < 0{e~^^^\loge\~^^'^). The Neumann condition at the inner 
boundary of ^ is a direct consequence of the assumption / e H^{A). □ 

Since the asymptotic behavior of the energy E^^ is the same as that of E'^^ (see (|2.2p ) for any 
|a;| < 0{e~^^^\ logej"'^'''''), most of the estimates proven for the density profile g hold true for g^^ as well, 
provided the phase ui satisfies the estimate required in Proposition |3?T] We sum up such estimates in the 
following 

Proposition 4.2 (Estimates for g^). 

Ifn^e-^\logs\-^ and IcjI < 0(e-5/4|iQg^|-3/4) 



12 TF I 



<0(e-V4|ioge|-3/4), (4.9) 

gliRrr.) - ll.9.|lio.(^) < ||p^nL==(B) {l + Ois'/'lhgel-'/^)) . (4.10) 
Moreover for any r £ A such that Ryi + 0{e\ logej"^) < r < max[i?,n, 1 — £'^/^| logep] 

|£M -P^^MI < 0{s-'^^\loge\-'^^) < 0(eV4|iog£|V4)^TF(^) ^ ^tf^^)^ (4^1^) 

and the maximum position i?ni('^) o/ 3^ satisfies the bounds 

Rn.{u;) > 1 - 0(£«/«| \oge\y'), \\9jUB\BiR„.)) < log^r/'). (4.12) 

Finally for any r such that r<Rh- 0{e'^/^), 

gl{r)<Ce-'\\oge\-' exp{--^}. (4.13) 

Proof. The results are proven exactly as the analogous statements contained in Propositions 12. 3[ 12.41 12.51 

andES □ 



4.2 Energy Decoupling and Optimal Phases 

The first step in the proof of the absence of vortices is a restriction of the GP energy to a subdomain of 
B and its splitting in a suitable energy functional plus the giant vortex profile energy. More precisely we 
consider the annulus A defined in (j4.ip with an inner radius i?< = i?h — e^^^ suitably chosen in such a 
way that outside A the estimates p.33p and (|4.13p yield the exponential smallness in e of both 'f'^^ and 
the density profile g^^. 
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We also recall the functional f introduced in (|4.4p . which is going to give the energy of the giant 
vortex profile, and the reduced energy 

£4v] := ^df 52 ||Vt;|2 - 2B^ ■ (tvyv) + e-'gl{l - \v\^f] , (4.14) 

where 

{iv,\/v) -.^ ^i{v\/v* - v*Vv) . (4.15) 
Proposition 4.3 (Reduction to an annulus). 

For any a; G Z such that \uj\ < 0{e~^^^\ loge\~^^^) and for e sufficiently small 

El^ + E^[u^] ~ 0(e°°) < < El^ + 0(e°°), (4.16) 

where the function is defined in A by the decomposition 

*'^''(^) ^■■gu{r)u^ir^e^p{tm-u)^}. (4.17) 

Proof. As in |CRY[ Proposition 5.4] the only ingredients for the proof of the above result are the ex- 
ponential smallness p.33p of outside A and the variational equation solved by g^. Note that the 
function is well defined away from the boundary dB where both v]/GP g^Yid g^ vanish. □ 

The idea behind the decomposition (|4.17p is that, if the phase factor w is chosen in a suitable way, 
the function u^^ obtained by the extraction from vJ/GP of a density g and the giant vortex phase, i.e., the 
phase factor exp{i([ri] — contains basically no more vorticity and \u^\ ^ 1 in some region close to 

the boundary of the trap. The optimal giant vortex phase is determined by inspecting the dependence 
on u} of the energy -E^^, i.e., one needs to identify the luq minimizing E^^ . 

Proposition 4.4 (Properties of the optimal phase ujq and density guig)- 

For every e > there exists an G Z minimizing E^ . Moreover one has 

0.0 = ^(1 + 0(|logsr^), J^drgl - j ^ o{l). (4.18) 

Proof. The existence of a minimizing ojq E Z, can be deduced as in |CRY[ Proposition 3.2] as well as the 
second estimate in (j4.18p . 

The estimate of luq is a straightforward consequence of the estimate^ on g^^ contained in Proposition 
14.21 since one has (recall the definition of the annulus ^buik in (|1.30p ) 

n I df(r-2_i)£^<t7(l + 0(£i/4|iog,|7/4)\ j df(r-2-l)pTF+o(,-i) f df55^< 

n(l + 0{e^'^\\oge\'''^)] f (^{r-^ -l)p^^ +0{e-'\loge\-^), (4.19) 

where we have used the fact that \A \ .4buik| < ^(£1 loge|^^) and the estimates (|4.10p and (|4.1ip . which 
also imply that 

sup g^^ir) <0{e-'\loger^). 

reA\Abuik 

On the other hand since 

n f df - 1) ^ ^ [1 _ ^4 ^ 2Rl\ogR^'] ^^{1 + 0{e\ \oge\)), (4.20) 

*Note that the second estimate in 114.1811 allows to extract the simple bound \lU()\ < 0{e~^) which guarantees that all 
the estimates proven in Section 14.11 apply to gujQ. 
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and 

/ dfr-'gl, > R-^^{1 - 0{e^^')) > 1 - 0(e| log£|), 
J A 

the result easily follows. □ 
The analogue in the whole ball B is discussed in the following 

Proposition 4.5 (Optimal phase Wopt)- 

For every e > there exists an Wopt G N fulfilling 

a;opt = :r^(l + 0(|log£r^) (4.21) 

which minimizes , i.e., 

E<^^^Egl^^. (4.22) 

Proof. The existence of Wopt can be proven as in Proposition l4.4l above. Moreover, as in }CRY[ Proposition 
3.2], it is not difficult to show that the following estimate 



dfg%, ( n ) = 0(1), (4.23) 



holds true, where ^opt is the minimizing density associated with Wopt- 

In order to extract the same information as in the proof of Proposition l4.4l one needs however to restrict the 
above integration to a domain comparable to A and this requires some further analysis of the properties 

Offfopt- 

Using a regularization of g^^^ as a trial function for S^^^^ and exploiting the exponential smallness (|4.13p 
one can easily show that 

E^^ =ES:^^^<ES:,+0{en, (4.24) 

which guarantees that all the estimates proven in Proposition 14.21 apply also to 5opt- Hence one can 
use the exponential smallness of gopt (see (|4.13p ) to estimate the integral inside B \ A, but this is not 
completely sufficient because the potential -B^opt contains a singular term at the origin ~ r^^ and one 
needs an additional estimate showing that ^opt vanishes as r — >■ 0. This is proven in Lemma |4. II below. 
By using (|4.26p and the analogue of (|4.13p . one thus obtains from (|4.23p 

J/r9%.{n^B_!^^=Oill (4.25) 

which implies the result exactly atj^ in the proof of Proposition 14.41 □ 
Lemma 4.1 (Pointwise estimate of gopt close to the origin). 

The density ijopt minimizing the functional f 5Ipt defined in (|1.35p satisfies the pointwise estimate 

5opt(r)<||5optlL^(H)(2r)[^'/2]. (4.26) 

for any < r < 1/2. 



^Note that the other estimates of gopt (analogous to those stated in Proposition 14. 211 which are needed to complete the 
proof can be derived from the energy bound I I4.24II . 
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Proof. The function W{r) := |l,9opt|loo (2r)[^/^l is a supersolution in [0,1/2] for the variational equation 
solved by ^opt, i-e., 

-Apopt + {[^] - i^opt)^r~^gopt - 2ft{[n] - ujopt)gopt + 2e"^gopt = Atoptffopt, 

since 

-AW+ m - uj,pt)^r-^W - 2n{[n] - uj,pt)W + 2e-^W^ - ^loptW > 

{ [{[n] - ujopt f - - cn] r-^ - 2n[[n] - ujopt) - Mopt} w{r) > cnwir) > 0, 

where we have used the estimate ^opt = ~^^(l~o(l)) and the fact that we are in the interval r g [0, 1/2]. 
Since at the boundary dBi/2 one has gopt(l/2) < ||.9opt|loo = W{l/2), the maximum principle (see, e.g., 
p]]) guarantees that gopt{r) < W{r) and therefore the result. □ 



4.3 Estimates of the Reduced Energies 

The next crucial step in the proof of the absence of vortices is the lower bound for the reduced energy 
functional and in the rest of this section we will focus on such a problem. Since the optimal phase 
ujQ as well as the associated density g^^^ can be fixed throughout the rest of the proof, we simplify the 
notation for the sake of clarity and set 

- £[v], g, R„,{ujo) i?m, B^„{r) B{r) = [fir - {[il] - cjo) r'^] e^, (4.27) 

and 

:=^drV{|VH' + e-V(l-|H')'}, (4-28) 
where i?< := i?h - e^^'^ and (see (|i?T|) l 

A := {fe B: r> i?<}. 
We also recall that u := u^^ is defined inside A by 

«'GP(f) =: g{r)u[r) exp {i([f}] - lo^)^} , 
and the annulus .4buik is (see (|1.30|) ) 

Auik ■.= [r(zB: i?> < r < 1 - e^'^\ logepj 

with i?> := i?h + e| loge|~^ (see (|1.3ip ). Note that thanks to the pointwise estimate (j4.1ip . we have the 
lower bound 

g^{r) > -f, ^ on Auik. (4.29) 

£| logel'^ 

We can now state the main result in this section, which is going to be the crucial ingredient in the 
proof of the absence of vortices: 

Proposition 4.6 (Bounds on the reduced energies). 

// O = r^oe^^l loge|^^ with > 2(37r)^^, then for e small enough 
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The proof of the above results is quite involved and before the discussion of its details, which is 
postponed to Section [1751 we are going to give a quick sketch of it together with the statement of several 
preliminary results. 

The main trick in the estimate of the reduced energy is an integration by parts of the second term 
in (|4.14p . which is made possible by the introduction of a potential function F{r) already considered in 
[CRYj . Such a function satisfies the key properties 

V-^F = 2g^B, F{R<) = 0, (4.31) 

and it is explicitly given by 

F{r) -.^2 J ds g^(s) ^Sls - ([SI] -loq)^ = 2 j dsg^{s)B{s) ■ e^,. (4.32) 

Other important properties of F are formulated in the next lemma and are basically straightforward 
consequences of (|4.18|) and the bound 

\B{r)\ <0{e-^) on A, (4.33) 

which follows from the definition of A. 

Lemma 4.2 (Useful properties of F). 

Let F be defined in (j4.32p . The following bounds hold true: 

WFh^iA) < ^^(e^')' \\^F\\L-iA) < (^i^^'l loger')- (4-34) 

Moreover one has the pointwise estimates 



-i?<|g2(r), z/re[i?<,^J, 
l + e-'\l-r\g^{r), z/r€[i?„„l]. 



|F(1)|<0(1), |F(r)|<C<i: ^'VL. . (4-35) 



Proof. Most of the proof follows from [CRYl Lemma 4.1]. The estimate (|LT8l) yields |F(1)| < 0(1). The 
last inequality in (|4.35p for r E [Rm, 1] is a consequence of this bound together with the identity 



1 

2/ 



F(r) = F(l)-2 / dsg'is)Bis)-e^, 

J r 

and the fact that g{r) is decreasing for r e [Rm, !]■ D 

Due to the lack of control of the behavior of the function u at the boundary dB, we need to use a 
suitable decomposition of F: An integration by parts (Stokes theorem) of the second term in (|4.14p would 
indeed give 



- 2 / dfg^B ■ ( m, Vw) = / dr _F(r)curl(m, Vu) — / da F{l){iu,dT-u), (4.36) 
J A ' J A JdB 

and the last term in the expression above clearly depends on u at the boundary. While Neumann boundary 
conditions allow to extract some information about u on dB and in particular an upper estimate for that 
term, on the opposite, if Dirichlct conditions are imposed, u is not even well posed on dB, since both 
^"^^ and g vanish there. A way out to avoid such a problem is the decomposition of F into a function 
vanishing on dB and another one whose gradient can be explicitly controlled: More precisely we set 



-Fout(r-) F(l) 



dss-V(s)l r dss-'g^is), (4.37) 
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so that 

V (.g^'VFout) = 0, -Fout(l) = (4.38) 

If we now define 

Fi„(r) :=i^(r)-Fout(r), (4.39) 

one can easily verify that 

V (5-'VFi„) = 2V • B{r)er, i^i„(l) = 0, (4.40) 
and, integrating by parts only the term involving Fi^ in (j4.36p we obtain 

-2/ di^g'^B ■ {iu,Vu) = - d?^ V-^i^out ■ («w, Vw) + / dr ^i„(r)curl(m, Vu). (4.41) 
Ja J a J a 

The energy £[u] can thus be rewritten as 

E[u]= j df |Vu|VFi„(r)curl(m,Vu) +£-2/ (1 - |up)^| - y df V^Fout • (m, Vm). (4.42) 

The first three terms above are the most important ones and their estimate is the key result in the proof 
of the absence of vortices. The last term on the other hand can be estimate separately and one can show 
that it yields only a smaller order correction. 

More precisely the first two terms can be estimated in terms of the vorticity of u: As in |CRYj . if we 
suppose that |ti| ^ 1 except in some balls {i3(aj, t)}^.^ J C N, whose radius t is much smaller than the 
width of .Abuik, and we denote by dj the degree of u around Oj , 

/ dr _Fin(r)curl(m, Vw) ~ V27rFi„(aj)dj. (4.43) 
and, optimizing w.r.t. the radius t. 



dr.g^ |V^.|^ > ^2V(«,)M,|log ( ) ^Y.^g\a,)\d,\\ \oge\. (4.44) 



|2 

Hence 



£[u] gJ2^7:\d,\ (i5'(%)|log£| +i^in(a,)) , (4.45) 

and, if f^o > 2(37r)~^, the sum between parenthesis is positive for any dj in the bulk (see Section fA. 2 [) . 
which means that vortices become energetically unfavorable. Note that there is an important difference 
with the analysis contained in jCRYj since F is replaced in the expression above by Fin- This is basically 
the main effect of Dirichlet boundary conditions. 

The starting point of the reduced energy estimate is given by the following preliminary upper bounds: 

Lemma 4.3 (Preliminary energy bounds). 

//rj ~ e-2|loge|-i ase-)-0, 

-FM < 0{e-^), £[u] > -0{e-^). (4.46) 

Proof. See [CRYl Lemma 4.2]. □ 
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4.4 Vortex Ball Construction and Jacobian Estimate 

In order to construct families of balls containing all the vortices of u, we need to exploit some local energy 
bound on J^[u\. However the bounds (j4.46p are not sufficient for our purposes, since they imply that the 
area of the set where u can possibly vanish is of order logep, whereas the vortex balls method requires 
to cover it by balls whose radii are much smaller than the width of A, which is 0{e\ logej). 
As in |CRY| there is a way out to this obstruction in the localization of the energy bound (|4.46p . given 
by the decomposition of the domain into suitable good and bad cells: 

Definition 4.1 (Good and bad cells). 

We decompose A into almost rectangular cells An, n Cz N, of side length 0{s\ loge|), given by 

An:={reA: 7? G [n9, {n + 1)9[} , (4.47) 

where 9 2'k/N and N ^ £^^\ logej^^ is the total number of cells. Let < a < ^ be a parameter to be 
fixed later on. 

We say that An is an a-good cell if 

drg^[\Vuf+e-'g^l-\u\^y}<s-'-^\\oge\, (4.48) 

whereas inside a-bad cells the (strict) inequality is reversed. We denote by and the numbers of 
a-good and bad cells and by GSa and BSa the sets covered by good and bad cells respectively. 

By definition of bad cells, one has that (|4.46p immediately implies 

< e^+"\loge\-^T[u] < Ce-'^+"\\oge\-^ < N, (4.49) 

i.e., there are very few a-bad cells. Note also that the final estimate (|4.30p implies that there arc actually 
no bad cells at all. 

We can now construct the vortex balls inside good cells but, since the density has to be large enough, 
we need to restrict the analysis to the subdomain w4buik C A (see p. 301) for its definition): 

Proposition 4.7 (Vortex ball construction inside good cells). 

For any < a < ^ and e small enough, there exists a finite collection := {B{d.i, ft)}jgj of disjoint 

balls with centers Si and radii Qi such that 

1. {feG^„nAuik: > |loge|-i} cUe/S,, 

2. for any a-good cell An, Y.i,B,nA^^0 Bi = e|loge|~^. 

Setting di :— deg{w, dBi\, if Bi C .Abuik H GSa, and di = otherwise, we have the lower bounds 

Arg" \Vu\^ > 2^ - «) \d^\a\a,) \\oge\ {\ c"^^^) ■ (4-50) 

Proo/. See [CMl Proposition 4.2]. □ 

Given a suitable family of disjoint balls as in the above proposition, one can prove that in the a-good 
set the vorticity measure of u will be close to a sum of Dirac masses, i.e., 

curl(m, Vu) ~ 2tt di5 {f ~ Si), 
where 6(f— Si) stands for the Dirac delta centered at Si. 
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Proposition 4.8 (Jacobian estimate). 

Let < a < i and cj) be any piecewise-C^ test function with compact support supp((/)) C .4buik H GSa ■ 
Let also {Bi}^^j := {B{di, gi)}^^j be a disjoint collection of balls as in Proposition \4.7\ 
Then setting di := deg{u, dBi}, if Bi C .4buik H GSa, o,nd di = otherwise, one has 



2Trdi(j>{di) - I 



df (p curl(zM, Vu) 



GSc n^buii 



<C||V0||^^(es^)£^|log£r^J-M 



Proof. See [CRYl Proposition 4.3] 



(4.51) 
□ 



4.5 Completion of the Proofs 

The main goal in this section is the proof of Proposition 1321 which will lead to the proof of Theorem II .41 
As anticipated before, the first important step is an integration by parts of the second term in (|4.14[) . 
but, since it has to be performed cell by cell, it generates boundary terms living on the frontiers between 
good and bad cells. Such terms are artificial, since the cell decomposition has no physical meaning, and 
we want to avoid having to estimate them. 

As in |CRY| we introduce an azimuthal partition of unity to get rid of these terms (see also |CRYi 
Definition 4.2 and Eq. (4.69)]): We define a pleasant set PSa as the set generated by good cells such 
that their neighbor cells are both good (pleasant cells) , whereas the average set ASa is made of good cells 
with exactly one good cell as neighbor (average cells) . Finally the unpleasant set U PSa contains all the 
remaining good and bad cells (unpleasant cells). Denoting by N^, N^^, N^^ the number of pleasant, 
average and unpleasant cells respectively, it is not difficult to see that 

N^^ < ^N^ < iV, iV^ < 2N^ < N. (4.52) 

The partition of unity is given by two functions Xin(^?) and XoutW such that Xin(^?) + XoutW = 1 for 
any d G [0, 27r] and 

Ji, ii^eUPSa, Ji, ii^ePSa, 

^-'^^^^^\o, ii^ePSa, ^'"^'"^^^jo, if^E^P^., ^'-''^ 

Since both functions vary from 1 to inside an average cell, one can always impose the bounds 

IVxoutI <0(e-'|logeri), l^Xin] < O{e-'\loge\-'), (4.54) 

because the side length of a cell is cx e| loge|. 

In order to apply the jacobian estimate proven in Proposition 14.81 to the function (f> ~ XmFim whose 
support is not contained in .4buik but only in A, we also need a radial partition of unity: We define two 
radii as (recall that i?> — Rh + £\ loge|~^ as in (|1.3ip ') 

Rt^, ■.= l~e\loge\-\ i?-,, :=i?>+e|loge|-\ (4.55) 

and two positive functions ^in{r) and ^out(?') satisfying £,in{r) + £,out{r) = 1 for any f £ A and (recall 
), i.e., R^=R^- eV6) 

. ,^ fl, ifr€[i?<,i?>], or[l-e3/2|logep,l] 

[0, if r e [i?^^t,i?+t], 

^'"^ ^ ■ jo, ifrG[i?<,i?>], or[l-£3/2|iog,|2^1]_ 14.50 
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Thanks to (|4.55p we can also assume 

iVCoutl < 0{e-'\\oge\), IVCinI < 0{e-'\\oge\ 

We are now ready to prove the bound on the reduced energies: 



(4.58) 



Proof of Proposition \4-(>] 

For the sake of simphcity we denote by {Bi},^^j := {B{ai, gi)}^^j a collection of disjoint balls as in 
Proposition |4?7l whereas the subset J C I identifies balls such that dj 7^ 0. 
The starting point is an integration by parts as in (j4.42p . i.e., 

£[u]= J dr |Vu|VFi„(r)curl(m,Vit) +£^^9'' (1 - lup)^} - y" df V^Fout • (m, Vm). (4.59) 

The last term in the expression above is the easiest to bound: By using the explicit expression of 
i^ont, one obtains 



(iu, Vu] 



<cj drg^{r)\u\\Vu\<c(^sJ dfg'^\u\'^+S-^J d?V |Vw|^^ < 

C{S + 6-^T[u]) < CT[u]^/^, (4.60) 



where we have introduced a parameter S and chosen S = T[uY^^ (recall that J-[u] > 0). 

The remaining term in ()4.59p can be estimated exactly as in |CRY[ Proof of Proposition 4.1], with 
only one difference due to the presence of Fi^ instead of F: Since by definition the former vanishes on 
dB, we can get rid of all the boundary terms (see, e.g., jCRYi Eq. (4.86)]) and the final result is, for 
some parameters 7, S that we fix below. 



J df |Vit|^ + Fin curl(m, Vu)| > 



^6n(aj)Mi 



(1-7)U-" .9n«j)|loge| 1-C 



,l 0g|l0g£| 
|log£| 



l^in(«,)l 



(I-7) / dfUts'lVup- / drUt\Fin{r)\\S/u\^ + h-6) f drg^\Wu\^ 
JA JA J A 

drg^\u\^ -C\\oge\-^F[u]. (4.61) 



_ SL 

We can now choose the parameters a, 8 and 7 as follows: 

log I log £ I 



7 = 2(5:= 



log£| 



-log I log £| 

a = a- 



|log£| 



(4.62) 



where 5 is a large enough constant (see below). 

Using the properties of the function Fl(r) := |log£| — |Fin| proven in Proposition I A. 1 [ we have for 
f^o > (37r)^i 



1 „2 



.g^(a,)|log£|-|^^in(a,)|>C£"11og£r 



for any € y^buik, so that 

(1 - 7) - «) g\a,) |log£| {\ Ci^^i^) - |Fi„(a,)| > 

\g^{a,) |log£| - |Fi„(a,)| - ^^^(a^.) log |log £| > C£-i|log£r^ ( 1 - 



Clog [log £| 
|l0g£| 



> 0, (4.63) 
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where we have used (|4.29p . 

On the other hand for any r € supp(^out) either |r — i?<| < Ce| loge|^^ or |r — 1| < Ce\ loge|^^, which 
by the bounds (|4.35p imply that in the first case 

|Fin(r)| < C (I loge\~'g\r) + \F,^t{r)\) < C {\ \oge\-'g\r) + l) < C\ loger'g^r), (4.64) 

thanks to (j4.29[) . whereas in the second case 

1 



li^inWI < |F(1) - FoutWI +2 ds \B{s)\g'is) < 

J r 

\F{1)\( f dss-\/{s)] f dss-'g\s) + C\loge\-'g^ir)<C\\oge\-'g\r). (4.65) 



\Jfl,< / Jr 

Note that in this second case there is no need to assume that r > i?,„ in order to use that g is decreasing 
in [r, 1]: By the bounds (j4.10p and (|4.11l) . for any 1 - e| loge|"^ < r <1, 

9\r) > (1 - o(l))pTF(i _ e\ log£|-i) > (1 - o(l))g2(7?„,), 

so that we can always bound in the integrals g'^{s) by (1 + o(l))g^(r). 
In conclusion 

I^^in(r)|«g2(r), (4.66) 

for any r S supp(^out) and thus 



drUt [(1 - 7) 9\r) ~ \F,,,{r)\\ \Vu\^ > 0. (4.67) 
Finally wc have from (gUT)), and 

df iVuf + Fi„(r) curl(iM, Vw)| > 

ci^|H^/d.V|V.P-C^ii^^/ d.V|.p^C|loge|-^M, (4.68) 



and adding 



dffj (1-1.0^ 

A = 



to both sides of (|4.68p and using (|4.59p and (|4.60p . we get the lower bound 

I llog^l enog\\oge\JuPs^uAS^ ^iij' 

valid for e small enough and JIq > (Stt)^^. But .g^|up = |\I/GP|2 < (j£^^\ loge|^^, whereas the side length 
of a cell is 0{e\ loge|), thus 

r ^ C\UPS U AS^l ^ ^^^^1 ^^^uP ^ ^A^ < Ce'+^T[u], (4.70) 

JuPScUASc e|ioge| 

by and Therefore becomes 

cr 1 ^ f l0g|l0g£| ^ Uogg| .aT-r 1 x-r ll/2\ (A7^\ 

I loge log log e L J L J J 
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Recalling the choice of a in (|4.62p . we now take a constant a > 2, so that 

|loge| ^ lloger-" ^ log|log£| 



log I log £| log I log £| |loge| 
and 

0{sn > e[u] > (^^l^-^N - ^uV^') (4.72) 

which yields both results. □ 



The proof of the energy asymptotics is essentially a corollary of the reduced energy estimates together 
with the discussion contained in Section 14.11 

Proof of Theorem \1.4\ 
By (ril) and KM 

> +£[u] - 0{en > Ell - C{\og I loge|)-2| \oge\\ 

but one can easily show that > E^^ — 0{e°°) by simply testing the functional on a suitable 
regularization of gi^^ and thus 

> ESI - C-Qog I loger^l log£|2 > ijgv _ c-Qog | log£|)-2| logep, 

which concludes the lower bound proof. 

The upper bound p.38|) is trivially obtained by testing the GP functional on a giant vortex function 
with phase [ft] — utopt (sec Proposition 14. 5p . □ 

Using the equations satisfied by vpG^P q^yi^ one can derive an equation satisfied by u: 
~V{g^Wu) - 2ig^B + - l) w Xg'^u 

where A = — /ig,^ . A useful estimate on the gradient of u follows from this equation and allows to 
conclude the proof of Theorem 11.31 We state the estimate for convenience and refer to jCRYi Lemma 
5.1] for its proof. 

Lemma 4.4 (Estimate for the gradient of u). 

Recall the definition of u in (j4.17p . There is a finite constant C such that 

l|V"ll.==(.A.„.)^^^7^- (4.73) 

We now complete the 
Proof of Theorem \1.3[ Suppose that at some point tq € .4buik we have 

||"(ro)|-l| >e'/'*|logep. 
Then, using (|4.73p . there is a constant C such that, for any re B{fo, Ce' loge^/^), we have 

\\u{r)\~\\>\e^^\\oge\\ 

This implies (recall (|4.29p ) 

dr4(l-|"^'>C|log£|^ 
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and thus 

T[u]>C\\ogs\^ (4.74) 

which is a contradiction with (|4.6D . 
We have thus proven that 

||*°^P-g^|<5^|H^-l|<cfe^ (4.75) 
on y^buik- The resuh then fohows by combining (|4.1ip and (|4.75p . □ 



Theorem 11.51 fohows as a corohary: 



Proof of Theorem \1.5l 

Given any i?> < r < 1 — £'^/^|log£p, Theorem 11.31 guarantees that deg{u,dBr} is weU defined and 
independent of r. Moreover one has 

2TT\deg{u,dBr}\< [ ds\u\-^\dru\<C [ ds\dru\, 

because u is bounded below in „4buik as a consequence of the proof of Theorem 11.31 Now integrating in 
r from i?> to 1 — e^/'^ \ logep both sides of the above expression and using the fact that the degree is 
independent of r because u has no vortices, we obtain 

27r|deg^,9SJ| < Ce-^\\oge\~^ f dr \Vu\ < Cs-'\\ogs\-' lAuikl'/' W^^WmA..,,) - (4-76) 

where we have used Cauchy-Schwarz incquahty and the fact that |w4buik| = 27r(l — £'^/^| logep — i?>) = 
0{e\ loge|). On the other hand, and imply 

\\^^\\mA....)<Cs'^'nogs\'^'- 

We conclude 

27T\deg{u,dBr}\ < C\logef 
and final result is thus a simple consequence of the definition (|4.17p . □ 



5 Rotational Symmetry Breaking 

We first introduce some notation that will be used in the proof of Theorem 11.61 The result stated there 
is equivalent to prove that no GP minimizer is a symmetric vortex, i.e., a wave function of the form 
/(r) explmi^}, n € Z. We therefore denote by En the energy obtained by minimizing the GP functional 
on symmetric vortices, i.e., 

E„:= inf £'^P[f{r)cxp{ind}]^£'^^[fn{r)c^p{ind}], (5.1) 
where /n(f) is the unique real minimizer. 

We also define n G N through miiini£zE„ =: Eni Note that a minimizing n certainly exists for any e 
thanks to the convexity in n of the functional. However such a minimizer needs not be unique because 
of some accidental degeneracy (there are at most 2 minimizers), which can be removed by a infinitesimal 
change of e. 

The next lemma contains several useful properties of fn: 
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Lemma 5.1 (Symmetric vortex minimizer). 

For any £ > and O 3> e^^, there exists some n G Z minimizing En and it satisfies the estimate 

n = n{i + o{e-^n-^)). 

The associated minimizer fn {r) is unique and, up to multiplication by a constant phase factor, it is given 
by a positive radial function vanishing only at r = and r = 1. Moreover it has a unique maximum at 
some point < i?* < 1 and satisfies the estimate ||/fi||^2(g\B^ j = o(l). 

Proof. Wc first notice that by setting fi —: [fi] — w for some cj S Z, one can easily recover the coupled 
minimization problem studied in Proposition 14.51 (see also Proposition 14. 4p for some different angular 
velocity fl. It is very easy to realize that the existence of a minimizing a; (and thus n) as well as the 
estimate w = 0{e~^) can be deduced in the same way as in Proposition 14.51 

On the other hand for any given n € Z the uniqueness and positivity of the minimizer /,i(r) can be 
deduced by standard arguments, whereas the existence of a unique maximum at some point < i?* < 1 
can be proven by a rearrangement argument as in Proposition 12.21 bv noticing that the potential n^r~^ 
is strictly decreasing. 

In order to prove the estimate, we first notice that the fact that n = $7(1 — 0{e^^il^^) implies the 
upper bound 

E^KE^"" +Oie-^)+Oie'/^n^^''), (5.2) 

which is a consequence of the pointwise estimates (j2.5p and (|4.26p together with the bound |i?[o]-n('')| < 
0(£^^) + 0([fi] - n) for any f e A. As in (|2.19p the above estimate yields 

\fi^ - < 0{e~'/-'n'^') + 0(e-i/4l]5/4), (5 3) 

We can thus repeat the proof of the pointwise estimate (|2.2ip and the final result is 

\g^{r) - pTF(,.)| < 0(eV2f^i/2) + o{e'/^n'^^), 

for any + e~^r2~^| loge|~^ < r'^ < I — loge^/^. The argument described in Remark [3] 
therefore gives 

Rl>l-o{e-'n-'), (5.4) 

which in addition to /? < 0{eil) implies the result. □ 

The main tool in the proof of the breaking of the rotational symmetry is the investigation of the 
second variation of the GP energy functional evaluated at some local minimizer: Given some 5* solving 
the variational equation 

- A«'-20-L* + 2£-2|*|^* = (5.5) 

where fiq, := £^^[^] +e"^||*l||, and some perturbation S(r) € H^iB), one has £^^[^! + eS] = £'^^[^] + 
e^Q*[S] +0{€^), where 

Q*[S]:= J df{|VS|'-2S*f1-LS + 4£-2|^ns|2-/x*|S|2}+2e-23fi J dr {^ii*)^E\ (5.6) 

By definition, if there exists some S e Hq{B) such that Qa([S] < 0, the associated local minimizer ^ is 
globally unstable and in particular can not be a global minimizer of the GP functional. 

Proof of Theorem \1.6[ 

Assuming that the GP minimizer was given by a symmetric vortex fn{r) exp{mi9} for some n, we explicitly 
exhibit a trial function S(r) such that the quadratic form Q*[S] evaluated at ^'(r) = fn{r) explihi)} is 
negative (for simplicity we denote it by Qn), which yields a contradiction with the assumption that the 
symmetric vortex is a global minimizer. 
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For any d > 1 we set 

S(f) {A{r) + B{r)) e*("+'^)'' + {A{r) - B{r)) e*(""^)*, (5.7) 

with 



0, if i?* < r < 1, ' I nRifnir), if i?, < r < 1. 



A trial function of this form was first introduced in |Seir[ Theorem 2] to prove symmetry breaking for a 
special class of trapping potential but here we replace in the original definition |Seir[ Eq. (2.30)] d with 
—d. Moreover in order to satisfy the condition S e Hq{B), we have modified the function A setting it 
equal to for r > R^,. Note that the function certainly belongs to H^{B) since fn is differentiable and 
A{r) + B{r) ^ r'^+<^ as r — )- 0, but B ^ Hq{B) because of the singularity in the derivative at r = i?*. 
We can simply borrow the explicit computations from |Seir[ Eqs. (2.31) and (2.33)] (recall that in our 
case there is no external potential, A vanishes for r > R^, d has to be replaced with —d and Q with 2Q) 
and, denoting by the chemical potential associated with fn, we obtain 

Q^[S] = 8^ dr r''+^Mr)f:-Ar) | (d + l)/ift - ^^^^/| + 2nn| + 47^2^^ dr ^fl{r). (5.9) 

Using the estimate (|5.3p . one immediately obtains that ^'^^ = = — fi^(l — o(l)). We can thus 
estimate the quantity between brackets in the first term in (|5.9p as 

{d + l)fin - 2e-^{d + 1)/| + 2nn < -fl^ {d-1- o(l)) (5.10) 

which implies the bound 

'''' dr r'''+'Mr)fUr) {{d + ~ 2e~\d + 1)/| + 20n} < 

- %^VL\d - 1 - 0(1)) dr r^'^+VfiM/ar) < 





12/ 



Rl^^^!l{R.) - {2d + 2)Rl'' f * dr r/|(r 



< 



Anfl^id - I- 0(1)) 



where we have used the fact that / is increasing between and i?, and the lower bound (|5.4p 
On the other hand the last term in (15.91) can be bounded as 



Ann'id - 1 - o{l))Rl''+^fi{R,) (1 - Cde-'n-^) , (5.11) 



Ann'd'Rl'' f dr fl{r) < 2n'd^Rl''-^ \\h\\l^iB\B^,) < oml^-'n^d' . (5.12) 

Hence (|ETT|) and ([5l^ yield 

Qf,[E] < -Airdn^Rl'^+'^fliR,) (l - d'^ - o(l) - Ce-^n'^d) < 
for any finite d> 2 and e small enough. □ 
Remark 5.1 (Flat Neumann case) 

The above proof applies with minor modifications to the case of the bounded trap B with Neumann 
conditions at the boundary dB: It is indeed sufficient to make the replacements in the trial function S 

B{r)^nr''Uir), 
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for any rGB and compute 

rR, 

Jo 

Now since fn is increasing in the Neumann case, /? < 0{eft) and fif, = —0^(1 — o(l)), the quadratic 
form can be made negative for d large enough. 



Appendix A 

In this Appendix we discuss some useful properties of the TF-like functionals involved in the analysis as 
well as the critical angular velocities. 

A.l The Thomas-Fermi Functionals 

The minimization of the TF functional introduced in (jl.l7p has already been discussed in other papers 
(see, e.g., }CY[ Appendix] or |CRY[ Appendix A]), so we only sum up here the main results: The 
minimizer among positive functions is unique and explicitly given by 

where [ • ]+ stands for the positive part and := E'^^ + e^^ ||p'^^||2- If ^ 2(-y7Fe)~^, the chemical 
potential is given by fi^^ = —il^Rf^ with 



and the TF minimizer can be rewritten as p^^{r) = ^e^Q^ [r^ — ^h]+' which makes explicit the fact 
that it vanishes for r < i?h- 

The corresponding ground state energy can be explicitly evaluated and is given by 

Note that above the second critical velocity, the annulus A^^ := supp (p^^) has a shrinking width of 
order e|loge| (see (jA.2p ') and the leading order term in the ground state energy asymptotics is —fl^, 
which is due to the convergence of p^^ to a distribution supported at the boundary of the trap. 
In the giant vortex regime another TF-like functional becomes more relevant, i.e., 

£r[p] ^ df {~n\'-p + Bl{r)p + e'^ p'] = df [{[il] - Lof r-\ + e-^'} - 2fi[17 - c.], (A.4) 

where the potential B^^ is defined in (|1.36p . uj E I1 and we have used the normalization in L^(B) of the 
density in the last term. The minimization of such a functional was studied in details in |CRY[ Appendix 
A] and we recall here only the most important fact, i.e., the ground state energy Ep-^ satisfies the estimate 



^TF ^ ^TF 



97re2 



O(e-2|logs|-i), (A.5) 



which suggests that it is minimized by a phase lo^^ 2(3\/7r£) ^ 
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A. 2 The Third Critical Angular Velocity Q 



C3 

In this last part of the Appendix we state the estimate of the critical velocity f2c3, which is defined 
as the angular velocity at which vortices disappear from the bulk of the condensate. To estimate this 
velocity we need to compare the vortex energy cost ^g'^g{r) \ loge| with the vortex energy gain |i^in(r)| (see 
(|4.39p . (|4.32p and (|4.37p ). In |CRY[ Appendix] a similar comparison is performed when the density g^^ is 
replaced by p^^ and it is shown that, if £7 > rt^s in the sense that f2 = 51o£~^| loge|^^ with > 2{3tt)~^, 
then the function 

where 

satisfies the lower bound 



F"(r):=2/ ds B^T.{r) ■ p'' (r), (A.7) 
H^^{r)>Ce-'\\oge\-^>0 (A.8) 



for any r such that r > i?> = i?h + £| loge] ^■ 
The analogous result for the original function 



Hir):=lgl{r)\\ogs\-\F,^{r)\, (A.9) 



is proven in the following 



Proposition A.l (Third critical velocity fica)- 

// f^o > 2(37r)^^ and e is small enough, there exists a finite constant C such that 

H{r) > Ce-'^\loge\-'^ > 
for any f such that r > _R> = i?h + £ \ logej^^. 

Proof. The result can be proven in the same way as |CRY1 Proposition A. 2] by noticing that |i^out| < 
|-P'(1)| < C'(l) and using such an estimate to replace Fin{r) with F{r) in H{r). □ 
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